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摘要 :【 目 的] 筛选 验证 意大利 蜜蜂 Apis mellifera ligustica Ж ё, DNA 3E28 6 5 di € JA 8 X 89 
SNP, [2 ik ] ЖР SERE E 0 Ascosphaera apis WPiA T A. 43 Je 8 E KR] З 日 龄 幼 
X ,根据 是 否 存 在 白垩 病症 状 进而 筛选 出 抗 病 个 体 和 易 感 个 体 。 基 于 前 期 重 测序 结果 中 意大利 蜜 
蜂 第 2 和 11 号 染色 体 DNA 非 编 区 与 抗 白垩 病 相 关 的 SNP 信息 ,利用 PCR 测序 的 方法 筛选 并 验证 
意大利 蜜蜂 幼虫 第 2 和 11 号 染色 体 DNA 非 编码 区 与 幼虫 抗 白垩 病 相 关 的 355 个 SNP。[【 结果】 发 
现 位 于 意大利 蜜蜂 第 11 号 染色 体 LOC100578413 基因 5 端的 非 编 码 区 的 SNP(T14570310C ) 在 抗 
病 个 体 中 了 等 位 基因 频率 高 于 C 等 位 基因 频率 , 且 抗 病 个 体 中 的 了 等 位 基因 频率 显著 高 于 易 感 幼 
下 中 的 了 等 位 基因 频率 ,表明 该 SNP 位 点 与 抗 白 垩 病 相 关 。 该 分 子 标记 对 抗 性 个 体 和 易 感 个 体 的 
判断 结果 与 前 期 第 选 的 编码 区 SNP( C2587245T) 分 子 标记 的 结果 一 致 。【 结 论 】 筛选 并 验证 意大利 
蜜蜂 第 11 号 染色 体 DNA 非 编 码 区 的 SNP(T14570310C ) 与 抗 白垩 病 相 关 。 该 位 点 为 抗 和 白垩 病 分 
子 辅助 选 育 提 供 新 的 分 子 标记 ,在 意大利 蜜蜂 白垩 病 早 期 检测 和 培育 白垩 病 抗 性 的 蜂 种 方面 具有 
关键 词 : 意大利 蜜蜂 ; 9 X38. 染色 体 DNA; 非 编 码 区 ; SNP; 抗 病 个 体 
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Screening and validation of SNPs associated with chalkbrood resistance in 


non-coding regions of DNA of chromosomes of Apis mellifera ligustica 

NIE Hong-Yi, XU Shu-Peng, WANG Xue-Yan, САО Yan, ZHU Ya-Nan, SU Song-Kun ' ( College of 
Bee Science, Fujian Agriculture and Forestry University, Fuzhou 350000, China) 

Abstract: [ Aim] This study aims to screen and verify the SNPs associated with chalkbrood resistance in 
the non-coding regions of DNA of chromosomes of Apis mellifera ligustica. [ Methods] Ascosphaera apis 
spores were inoculated to the artificially reared 3-day-old larvae of A. m. ligustica, and the disease- 
resistant individuals and susceptible larvae were screened based on whether they displayed typical 
symptoms. According to the SNP information associated with chalkbrood resistance in the non-coding 
regions of DNA of chromosomes 2 and 11 from our previous resequencing results, 55 SNPs associated with 
chalkbrood resistance located in the non-coding regions of DNA of chromosomes 2 and 11 of A. m. 
ligustica larva were screened and validated by sequencing of PCR products. [Results] The results 
revealed that the SNP T14570310C located in the 5' non-coding region of Ше L0C100578413 gene on 


chromosome 11 of A. m. ligustica showed a higher T allele frequency than the C allele frequency in the 
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disease-resistant individuals and a significantly higher T allele frequency in the disease-resistant 


individuals than in the susceptible larvae, suggesting that this SNP is associated with the resistance to 


chalkbrood disease. This molecular marker was consistent with the SNP C25872457T located in the coding 


region screened previously in respect of the distinction between the disease-resistant and susceptible 
individuals. [Conclusion] The SNP T14570310C in the non-coding region of DNA of chromosome 11 of 


А. m. ligustica was screened and verified to be associated with the resistance to chalkbrood disease. This 


SNP provides a new molecular marker for molecular-assisted breeding against chalkbrood disease, and has 


important significance in early detection of chalkbrood infection and breeding bees resistant to chalkbrood 


disease. 


Key words: Apis mellifera ligustica; chalkbrood disease; chromosomal DNA; non-coding region; SNP; 


disease-resistant individuals 





ЖИЕН EDS Н 2 ERRE A Ascosphaera. apis Jg 
染 蜜蜂 幼虫 所 引起 的 一 种 真菌 性 疾病 ,最 早 在 意 大 
ЖІЗЕШЕ Apis mellifera ligustica 中 发 现 ( Aronstein and 
Murray, 2010) , 且 对 意大利 密 蜂 蜂 群 群 势 和 生产 性 
能 影响 较 大 ( Reynaldi et al., 2015) 。 近 几 年 ,在 中 
华 蜜 蜂 Apis cerana cerana ( Chen et al., 2018) 、 木 蜂 
Xylocopa californica arizonensis ЖІ НЕ МЕ Bombus 


griseocollis 中 也 有 检测 到 白垩 病 ( Maxfield-Taylor et 


















































免疫 应 答 的 基因 和 通路 ,揭示 了 发 育 、 能 量 代谢 、 信 
号 转 导 ,疾病 与 免疫 应 答 相 关联 ,为 研究 蜜蜂 对 白垩 
病 免 疫 应 答 的 分 子 机 制 提供 了 有 价值 的 基因 转录 信 
息 (Xu et al., 2019) 。 

单 核 背 酸 多 态 性 (single nucleotide polymorphism, 
SNP) 由 于 其 数量 多 、 分 布 广 、. 易 于 自动 化 分 析 等 优 
势 ,被 认为 是 目前 应 用 前 景 最 好 的 遗传 标记 
(Invernizzi et al., 2011) 。Holloway 等 (2013 ) 对 不 同 
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抗 白垩 病 能 力 个 体 的 SNP. 位 点 与 抗 病 性 状 之 间 的 











(Heath, 1982) ,异性 菌 丝 通 过 相互 接触 产生 孢子 
( Anderson and Gibson, 1998; Piggeler，2001 ) 。 球 
均 菌 孢子 为 长 椭圆 形 (Li et al., 2018) ,通过 幼虫 进 
食 进入 到 肠 道中 ,在 围 食 膜 中 萌发 (Theantana and 
Chantawannakul , 2008) ,之 后 菌 丝 从 体腔 内 部 穿 出 
肠 壁 ,尾部 长 出 白色 阔 丝 ( 郊 志 阳 等 , 2011; 李江 红 
等 , 2012) ,死亡 后 继续 脱水 干燥 形成 白垩 状 干 己 
(Spiltoir, 1955) 。 由 于 白垩 病 发 病 早期 没有 明显 症 
状 , 一 旦 发 现 已 进入 晚期 (Flores et al., 20053) , H. 
孢子 存活 时 间 可 以 达到 15 4 (Flores et al., 
2005b)。 虽 然 蜜蜂 的 卫生 行为 可 以 将 发 病 或 死亡 
蜜蜂 拖 出 集 外 ,但 如 果 孢 子 已 经 形成 ,蜜蜂 的 卫生 行 
为 反而 会 促使 抱 子 传播 ,进而 加 剧 白 垩 病 病 情 
( Invernizzi et al., 2011) 。 白 垩 病 的 暴发 严重 危害 蜂 
群 群 势 , 对 养 蜂 人 造成 巨大 的 损失 ( 赵 红 霞 等 , 2014) 。 
随 着 蜜蜂 球 喜 菌 基因 组 测序 完成 (Qin et al., 
2006) , JIE JF. f XT 9E ER Е W HE Zhi fij. КМА (non- 
coding КМА, ncRNA) 的 研究 ,鉴定 了 379 个 长 链 非 
编码 RNA( long non-coding RNA, IncRNA) 和 118 个 
微小 RNA(microRNA, miRNA) ,不 仅 扩大 了 真菌 的 
IncRNA 数据 库 ,而 且 为 进一步 研究 潜在 的 IncRNA 
介 导 的 蜜蜂 球 圳 菌 基因 调控 奠定 了 基础 (Guo et al., 
2018). 。 通 过 转录 组 测序 技术 来 分 析 蜜 蜂 对 球 寺 菌 








































































































关系 进行 了 人 研究 ,鉴定 出 蜜蜂 幼虫 第 2 和 11 号 染色 
体 与 蜜蜂 幼虫 抗 白垩 病 性 状 显 著 相关 ,为 分 子 辅助 
选 育 提 供 了 基础 。 本 团队 受 励 民 等 (2014 ) 基于 重 
测序 技术 筛选 出 抗 白垩 病 蜜 蜂 幼 虫 特有 SNP. 位 点 
为 1 620 个 ,其 中 位 于 编码 区 与 蜜蜂 抗 白垩 病 相关 
的 SNP 位 点 680 个 (118 个 位 于 第 11 号 染色 体 上 ， 
52 个 位 于 第 2 号 染色 体 上 )。 刘 元 珍 等 在 前 期 研究 
的 基础 上 通过 实验 室 培养 蜜蜂 幼虫 筛选 抗 病 幼 虫 个 
体 和 易 感 病 幼 虫 个 体 得 选 验证 出 一 个 可 以 评估 蜂 群 
抗 白 亚 病 能 力 的 SNP 分 子 标记 C2587245T ,该 标记 
位 于 MRJPS 基因 第 2 个 内 含 子 区 域内 ,其 在 抗 白 亚 
病 蜂 群 中 C 等 位 基因 频率 (Pe) 显 著 大 于 易 感 病 蜂 
群 中 的 频率 , 且 C2587245T 位 点 Pe 高 的 蜂 群 白垩 病 
发 病 率 较 低 。 因 此 ,C2587245T 可 以 作为 一 个 分 子 标 
记 选 育 抗 白垩 病 的 蜂 群 或 蜂王 (Liu et al., 2016) 。 
除了 编码 区 序列 SNP 位 点 对 生物 性 状 表现 具 
有 重要 的 作用 外 ,生物 基因 组 DNA 非 编 码 区 也 存在 
与 疾病 相关 的 SNP 位 点 。 非 编码 DNA 可 以 转录 出 
非 编 码 КМА, ncRNA 可 通过 自身 转录 、DNA 甲 基 
化 组 蛋白 修饰 、 重 塑 染 色 质 等 途径 调控 基因 表达 ， 
影响 多 种 生物 学 过 程 ( 秦 丹 和 人 徐 存 挫 , 2013 ) 。 鉴 于 
DNA 非 编码 区 的 生物 学 意义 ,本 研究 在 前 期 已 有 的 
基础 上 ,筛选 并 验证 意大利 蜜蜂 幼虫 第 2 和 11 号 染 
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色 体 DNA 非 编 码 区 上 与 白垩 病 抗 病 相 关 的 SNP. 位 
点 ,为 白垩 病 的 诊断 提供 更 多 的 分 子 标 记 , 加 速 分 子 
辅助 筛选 和 培育 抗 白 亚 病 蜂 种 进程 ,具有 重要 的 实 
践 应 用 价值 。 








1 材料 与 方法 


1.1 实验 蜂 群 与 样本 采集 

本 实验 采用 的 5 个 实验 蜂 群 均 为 2017 年 来 自 
福建 福清 蜂 场 的 “ 蜂 强 1 号 ”意大利 蜜蜂 ,该 蜂 种 是 
采用 蜜蜂 闭锁 群体 育种 (closed-population breeding) 
方法 ,以 蜂王 浆 优 质 高 产 抗 白垩 病 为 目标 ,从 浙江 杭 
MERO OE ЯШ 慈溪 、 桐 廊 、 兰溪、 龙 游 . 江 山 、 
新 昌 和 德清 11 个 不 同 血统 的 蜂 场 ,选择 蜂王 浆 产 量 
高 .没有 白垩 病 发 病史 的 蜂 群 作为 育种 素材 ,分别 培 
育 雄 蜂 和 处 女王 ,采用 自然 隔离 交尾 和 蜜蜂 人 工 授 
精 技术 ,后 期 结合 分 子 辅助 选 育 技术 ,注重 抗 白 亚 病 
经 济 性 状 和 蜂王 浆 优 质 高 产 性 状 的 选 育 ,兼顾 蜂 群 
繁殖 速度 和 蜂蜜 产量 性 状 的 考察 选 育 ,经 过 多 代 的 
培育 ,形成 遗传 相对 稳定 、 抗 白垩 病 能 力 强 、 群 势 强 、 
蜂王 浆 优 质 高 产 的 蜂 种 。 实 验 蜂 群 饲养 于 福建 农林 
大 学 蜂 学 学 院 , 蜂 群 群 势 保持 一 致 。 为 了 获得 相同 
日 龄 的 幼虫 ,将 5 群 实验 蜂 群 的 蜂王 进行 限 王 产 卵 ， 
然后 从 每 个 蜂 群 分 别 取 192 头 2 日 龄 幼虫 转移 到 4 
个 48 孔 细 胞 培养 板 中 (预先 加 入 10 pL 幼虫 食物 ) 
共计 20 板 。 移 虫 过 程 中 ,每 移 取 10 头 幼虫 用 70% 
的 酒精 消毒 移 虫 针 ,防止 病原 交叉 传染 。 幼 虫 食物 
的 配方 ` 饲 喂 量 与 培养 方法 参考 前 人 报道 
( Crailsheim et al., 2013) 。3 НІР ІШІНЕ ЖЕ 
RAT, TIRY Vojvodie 等 (2011 ) 的 方法 培养 
和 纯化 蜜蜂 球 宫 菌 获 得 ,每 头 幼 虫 饲 喂 10 kL 含有 
SOME ER E TEL TL T YR BEI 5 x 10 个 /pL 的 食物 
(Jensen et al., 2013) 。 每 天 在 固定 时 间 饲 喂 幼 虫 ， 
从 6 日 龄 开始 收集 白垩 病 发 病 个 体 ( 发 病死 亡 特 
征 : 忠 体 不 蠕动 .身体 无 弹性 、 无 光泽 、 尾 部 有 白色 菌 
丝 ) 作 为 易 感 病 样本 ; 幼虫 到 达 10 日 龄 时 ,将 化 肾 
的 未 死亡 的 个 体 收 集 起 来 作为 抗 病 样 本 。 易 感 样本 
和 抗 病 样本 置 于 - 807C 冰箱 保存 。 
12 主要 仪器 与 试剂 

梯度 PCR 仪 (Applied Biosystems) , NanoDrop 2000 
《基因 有 限 公 司 )、 低 温 高 速 离心 机 (Eppendorf) 、 全 
自动 凝 胶 成 像 分 析 仪 (上 海 培 清 科技 有 限 公 司 )、 核 
酸 提取 仪 KINGFISHER FLEX ( Thermo Scientific) 等 。 
蜂王 浆 ( 购 自 杭 州 余杭 蜂 场 ,通过 低温 冷冻 运 

















































































































































































































输 , 且 储存 于 -80Y 冰箱 中 ) ;和 葡萄糖 (D-glucose)、 
果糖 (D-fructose) .酵母 提取 物 、 氧 化 钠 、 硫 酸 钾 、 蛋 
Fi sus ЕСЕТ ЕР TUER BE (3 43 RAST );2 x 
EasyTaq PCR SuperMix (北京 全 式 金 生物 科技 有 限 
公司 ); LabServ Cell & Tissue DNA Kit ( Thermo 
Scientific ) 等 。 
1.3 引物 设计 

根据 受 励 民 (2014) 的 蜜蜂 全 基因 组 重 测序 结 
果盘 选 得 到 了 关于 蜜蜂 幼虫 DNA 非 编码 区 与 抗 白 
垩 病 相 关 的 SNP 信息 ,其 中 位 于 第 2 和 11 号 染色 
体 DNA 非 编码 区 的 候选 SNP 分 别 为 24 个 和 31 个 
( 表 1)。 从 NCBI 数 据 库 中 搜索 并 下 载 这些 位 点 前 
后 各 500 bp 的 序列 ,利用 Primer Premier 6 设计 特异 
性 引物 ,PCR 产物 片段 大 小 约 为 400 bp ,包含 SNP 
位 点 , 且 正 反 向 引物 Tm 值 接近 。 
1.4 基因 组 DNA 提取 

利用 核酸 提取 仪 KINGFISHER FLEX, 并 按照 
LabServ Cell & Tissue DNA Kit 的 操作 步骤 提取 1. 1 
节 意 大 利 蜜蜂 幼虫 各 处理 样品 的 DNA。 主 要 步骤 
如 下 :用 灭 菌 后 的 眼科 手术 前 前 取 10 mg 实验 样本 
( 抗 病 与 易 感 病 样 本 分 开 操作 , 且 每 剪 1 头 后 ,用 
TOM 的 酒精 对 剪刀 进行 擦洗 ) , 置 于 1.5 mL 离心 管 
中 ,加 入 100 uL Lysis Buffer 和 5 uL 和 蛋白酶 涡 旋 
混 匀 ;55C 水 浴 2 hb, 期 间 颠 倒 混 匀 若干 次 以 提高 裂 
解 效 果 ,然后 室温 静 置 10 min ;按照 核 酸 提 取 仪 的 操 
作 说 明 运 行 DNA 程序 ;提取 DNA 的 浓度 和 纯度 利 
用 NanoDrop 2000 检测 。 
1.5 PCR 扩 增 和 测序 

PCR 扩 增 体系 (20 uL): 2xEasyTaf PCR 
SuperMix 10 pL， 上 下 游 引 物 (10 pmol/L) 各 0.4 
uL, 1.4 节 提 取 的 基因 组 DNA 2 uL, 加 ddH,O 至 
20 uL, PCR 扩 增 条 件 : 94C 预 变性 5 шіп; 94% 变 
性 30 s, 55*C3R Ж 30 в, 72% 延伸 1 min, 循 环 数 为 
35; 72% 再 延伸 10 min, 4% 保存 1 ho PCR 产物 由 
铂 尚 生物 技术 (上 海 ) 有 限 公司 进行 测序 。 
1.6 SNP 位 点 初步 验证 与 扩大 群体 验证 

随机 选取 意大利 蜜蜂 抗 病 个 体 和 易 感 个 体 各 
12 头 , 分 别提 取 其 ОМА 用 于 初步 验证 前 期 本 团队 
重 测序 获得 的 55 个 (2 号 染色 体 上 24 个 ,11 号 染色 
体 上 31 个 ) 与 抗 白垩 病 相关 的 DNA 非 编 码 区 SNP 
位 点 。 以 每 个 个 体 DNA 为 模板 ,利用 包含 SNP 位 
点 的 引物 ( 表 1) 分 别 进行 PCR( 同 1.5 节 ), 并 将 
PCR 产物 进行 测序 ,最 后 计算 每 个 SNP 位 点 等 位 基 
因 的 频率 。 我 们 将 初步 验证 中 抗 病 个 体 与 易 感 个 体 
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等 位 基因 频率 分 布 无 差异 显著 性 的 SNP 位 点 淘汰 ， 
而 将 等 位 基因 频率 分 布 存 在 差异 显著 性 的 SNP 位 
点 进一步 扩大 群体 验证 。 在 扩大 群体 验证 中 ,我 们 
再 随机 选取 意大利 蜜蜂 抗 病 个 体 和 易 感 个 体 各 90 
头 ,分 别提 取 其 DNA. 用 于 PCR 产物 测序 ,计算 每 个 
SNP 位 点 的 等 位 基因 频率 ,并 验证 这 些 SNP 位 点 的 
等 位 基因 频率 趋势 是 否 与 初步 验证 的 结果 一 致 。 同 
时 ,我 们 在 抗 病 个 体 和 易 感 个 体 中 也 检测 前 期 已 报 
道 的 编码 区 与 抗日 垩 病 相 关 的 SNP(C2587245T) 等 
位 基因 频率 ,验证 Pc 在 抗 性 个 体 中 的 频率 是 否 显 
著 高 于 其 在 易 感 个 体 中 的 频率 。 
1.7 生物 信息 学 分 析 

在 NCBI( https; // www. ncbi. nlm. nih. gov/) f 
酸 数 据 库 在 线 检索 SNP 对 应 的 NW 编号 ( 表 1 ) , 根 



























































据 NW 序列 与 染色 体 序列 的 对 应 关系 ,确定 该 SNP 
在 基因 组 上 对 应 的 染色 体 编号 及 位 置 ;NCBI ERIS 
索 染 色 体 的 编号 ,根据 SNP 在 染色 体 序 列 上 的 位 置 
定位 到 该 SNP 并 得 到 该 SNP 附近 的 基因 信息 。 
1.8 数据 分 析 

在 NCBI 意大利 蜜蜂 基因 组 数据 库 中 ,下 载 每 
个 SNP 位 点 附近 的 一 段 序列 作为 参考 序列 ,并 利用 
Lasergene (SeqMan) 软件 将 测序 结果 与 参考 序列 进 
行 校 对 分 析 , 排除 假 阳性 的 SNP 位 点 ,记录 下 SNP 
在 每 个 样本 的 基因 型 ( 纯 合 或 者 杂 合 ) ,统计 在 抗 病 
与 易 感 病 个 体 中 各 基因 型 的 数量 , 计算 基因 频率 。 
将 相关 数据 导入 IBM SPSS Statistics 21 统计 分 析 软 
件 ,利用 卡 方 检验 进行 差异 显著 性 分 析 。 


























R1 意大利 蜜蜂 幼虫 第 2 和 11 号 染色 体 ОМА 非 编 码 区 与 抗 和 白垩 病 相 关 的 候选 SNP 及 其 引物 序列 


Table 1 SNPs associated with chalkbrood resistance in non-coding regions of DNA of chromosomes 2 and 11 


of Apis mellifera ligustica larva and sequences of their primers 











gu SNP 所 在 染色 体 骨架 (染色 体 ) SNP 在 基因 组 上 的 位 置 
号 
Location of the SNP on the Position of SNPs in 

0. 


scaffold ( Chromosome ) the genome 


SNP 在 基因 组 上 的 序列 
Sequence of SNPs in 





引物 序列 (5' -37) 


Primer sequences 
the genome 








1 NW_003378123.1 (Chr. 2) 3 447 130 
9; NW 003378123.1 (Chr. 2) 1 916 903 
3 NW 003378123.1 (Chr. 2) 1 326 498 
4 NW 003378123.1 (Chr. 2) 11 528 
5 NW 003378111.1 (Chr. 2) 1 108 857 
6 NW 003378082. 1 (Chr. 2) 605 442 
7 NW 003378082.1 (Chr. 2) 605 427 
8 NW 003378082.1 (Chr. 2) 602 578 
9 NW 003378082.1 (Chr. 2) 265 918 
10 NW 003378082.1 (Chr. 2) 233 106 
11 NW 003378082.1 (Chr. 2) 92 706 
12 NW 003378058.1 (Chr. 2) 696 915 
13 NW 003378058.1 (Chr. 2) 144 846 
14 NW 003378058.1 (Chr. 2) 116 986 


: GCCGTCGCTTCTACTCTCTC 

: AGCACGACATTCCGTTGGAT 
: CCGCTATCTTTGTCGCCTCT 

: GGAATGTCTCCCCGTTCGTT 

: ССАСТССТТААСАСТАСАСА 
: СТААСТССТААСАССАСААСА 
: GACCAACCATTGCTTCTTC 

: GCTTGACCGTCCATAATCT 

: ССТТСАТТССТААТСАСАСТС 
: CTAATTGCGTTCGTGACAT 

: GACACAGAAGACCATTGAAC 
: САТСССТССТССААСААТ 

: GACACAGAAGACCATTGAAC 
: САТСССТССТССААСААТ 

: GAAAGAAAGACAGAAGAGAAGG 
: CTATCCTGACCTATCGTTACC 
: GATGATGACTTCGTGTGAAT 
: ATTAACCCTCTGCCGATG 

: ТТТАССАСССССССАААСАА 
: TACTGGAACAGGGGGATGGT 
: GCGGTGGTGTGATCTCTTGA 
: САТАААСССАСССССАССТС 
: TCTACTACGAATCACGGATG 
: GGAAGGTTGAAGGTAAGGAA 
: СТАСАТСТТССССССААССА 
: ТСТТССССАСАААССААССС 
: ACGCAATCCCTGTTCACCTC 

: GTAGCTCCTCACCGCAAAGT 


A 


G 


G 


G 
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续 表 1 Table 1 continued 








" SNP 所 在 染色 体 骨架 (染色 体 ) 
Location of the SNP on the 
No. 


scaffold ( Chromosome ) 


SNP 在 基因 组 上 的 位 置 
Position of SNPs in 





the genome 





SNP 在 基因 组 上 的 序列 
Sequence of SNPs in 


the genome 


引物 序列 (5 -37) 


Primer sequences 





15 NW 003378040.1 (Chr. 2) 


16 NW 003377991.1 (Chr. 2) 


- 


17 NW 003377991. 


(chr. 2) 


18 NW 003377991.1 (Chr. 2) 


19 NW 003377991.1 (Chr. 2) 


20 NW 003377991.1 (Chr. 2) 


21 NW 003377976.1 (Chr. 2) 


22 NW 003377965.1 (Chr. 2) 


23 NW 003377911.1 (Chr. 2) 


24 NW 003377911.1 (Chr. 2) 


25 NW 003378155.1 (Chr. 11) 


26 NW 003378155.1 (Chr. 11) 


27 NW. 003378155.1 (Chr. 11) 
28 NW 003378155.1 (Chr. 11) 
29 NW 003378155.1 (Chr. 11) 
30 NW 003378155.1 (Chr. 11) 
31 NW 003378088.1 (Chr. 11) 


32 NW 003378088.1 (Chr. 11) 


33 NW 003378088.1 (Chr. 11) 


34 NW 003378088.1 (Chr. 11) 


35 NW 003378088.1 (Chr. 11) 


36 NW 003378088.1 (Chr. 11) 





37 NW. 003378088. 


- 


( Chr. 11) 


569 303 


1 363 256 


1 293 878 


487 139 


3 844 


3 791 


403 466 


109 593 


6 538 


6 462 


980 174 


980 137 


980 111 


980 101 


160 793 


156 247 


4 009 533 


3 928 681 


3 156 753 


3 117 278 


3 068 820 


2 920 845 


2 597 336 


T 


чел Nou mou mom mom mom mom onm Got motu mon mom mom mom mon mou mouumw'usmumu's min mwimsmwi'mumw'em 


: AACAACGTCCGCCACTACTT 
: GCGGCAAGCAATAATACCCC 
: CCACGCAATCCTCTAACA 

: GCTTCGATCAAGGTCAGT 

: GGCTGGACCGCATAATCCAT 
: ТСАССААССССТССААСААА 
: ААСТТСТТССАСССАСТАТА 

: АСААТСТСАТТССААТСТСС 

: СССТТАТТССТССАСАСА 

: СТССАТССТСАТССАТАСАТ 
: СССТТАТТССТССАСАСА 

: СТССАТССТСАТССАТАСАТ 
: АСССТТАТСАСССТССТСТС 

: СААСССААТССССТСТСТТТ 

: ТССАТСТТТСАТССТССАА 

: TTCAGAGATCAGTGTGGAAT 
: TTGTCACTGTIGCTCGTT 

: СААТСТСАСТССААСААТСС 
: TTGTCACTGTTGCTCGTT 

: СААТСТСАСТССААСААТСС 
: GGTTATATCATTGTGGTCACTC 
: ССССАТССТТТСАТТСТТАТТС 
: GGTTATATCATTGTGGTCACTC 
: ССССАТССТТТСАТТСТТАТТС 
: GGTTATATCATTGTGGTCACTC 
: ССССАТССТТТСАТТСТТАТТС 
: GGTTATATCATTGTGGTCACTC 
: ССССАТССТТТСАТТСТТАТТС 
: СТАСАТТССТСССАСТТАТАС 
: GCGTITAGATCGTCTTGAAA 

: CAATCGCTATCTAACTTATCGG 
: ССТТСААСААССААСТАССА 
: АССАТТСТСАТТСТТССТТС 

: GCTCGGTTATTCAGTTGTC 

: СССССААСАТСТААСААТС 

: ТТАССАССАССССАСАТА 

: TTTCTTCTGGGATTCTGTGA 

: СТТССТТТСТТАССТАССАТС 

: ССАССАСАТААССААТАСАС 
: СССАССАСТТТСАТССАА 

: ТСАССТСТССАСТССАТС 

: ССААСААСССАААТСТТАСС 
: СТСТТСССТАСАСАССАА 

: АССААСАССАТААССТТСА 

: ААСАССТАССААССАСТС 

: СССАССААСТТСТТАССА 
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续 表 1 Table 1 continued 











Qu ЗАРРИНА СА) SNP 在 基因 组 上 的 位 置 
Location of the SNP on the Position of SNPs in 
o. 


scaffold ( Chromosome) the genome 





SNP 在 基因 组 上 的 序列 


01057 3 
Sequence of SNPs in 引物 序列 ( ) 


Primer sequences 
the genome 





38 NW 003378088.1 (Chr. 11) 2 594 038 
39 NW 003378088.1 (Chr. 11) 2 226 899 
40 NW. 003378088. 1 (Chr. 11) 1 880 352 
41 NW. 003378088.1 (Chr. 11) 1 440 682 
42 NW. 003378088. 1 (Chr. 11) 1 230 068 
43 NW. 003378088.1 (Chr. 11) 311 157 
44 NW. 003378054. 1 (Chr. 11) 1 284 280 
45 NW. 003378054.1 (Chr. 11) 1 060 121 
46 NW 003378043.1 (Chr. 11) 613 167 
47 NW 003377985.1 (Chr. 11) 13 310 
48 NW 003377973.1 (Chr. 11) 425 874 
49 NW 003377973.1 (Chr. 11) 425 835 
50 NW 003377973.1 (Chr. 11) 420 935 
51 NW 003377973.1 (Chr. 11) 323 967 
52 NW. 003377923.1 (Chr. 11) 476 412 
53 NW. 003377884. 1 (Chr. 11) 594 195 
54 NW. 003377884. 1 (Chr. 11) 497 716 
55 NW 003377857.1 (Chr. 11) 112 321 





2 结果 


2.1 意大利 蜜蜂 染色 体 DNA 非 编 码 区 SNP 位 点 

基于 前 期 本 团队 重 测序 获得 的 55 个 与 抗 白 亚 
病 相关 的 意大利 蜜蜂 第 2 和 11 号 染色 体 DNA 非 编 
码 区 SNP 位 点 ,我 们 通过 PCR 和 测序 方法 ,计算 每 





























: ТССААСААТССТСТСССТТТС 
: АТТССАТСТТСССАССССАА 
: СААСТСТСААСААСТТСААС 
: GGAGATTAGGAGCAGGAG 

: ССАСССТТССААТТААСТТ 

: СТТСТТСТТСССАСТАТАСА 


C 


: CCATGTAGAGCAGAGTGTT 
: GGTGATTCTCCGCAGTTA 

: СТТСТТСТТССТІТСТТССТ 
: ТСТТССССТСТТАСТСТТ 

: СТССТСТСАСТСАССТАТАС 
: TTCAATACCGTCGCTAGATT 


с 


: CCCCACACTATGCCACGTAA 

: САССАСТСТСТСТСТСТССС 

: ССАССАТААССТАТСАСААТАС 
: GCGACTAAGTAAGAGCTGTAT 
: АСТСАСАСААССАТААССАСАСТС 
: ТСАССТСТССАААССАТСТС 

: AAGATTCCGGTACTGTGCGA 

: GGACATTGCCTTGGGGACTT 

: CCAAGCATGACAGAGATGA 

: GCAAACATAAGTGGAGGATAC 
: CCAAGCATGACAGAGATGA 

: GCAAACATAAGTGGAGGATAC 
: СССАССТТССТТТСССТТТТ 

: GAATTGCTAGGCAACACGGC 


G 


G 


: CATTCTACTGCTTCACTAAC 

: GGTATCCTAATTGTTGATGG 

: CGACAGAGATAGGTATTCATC 

: САССТТАСАССТААСТАТАСТС 
: ACCCACGCACTTTCTCCTTT 

: AATCAACCGAACCTGCCGAT 

: TTACACGACGCCATCTAC 

: CATCTGTTGACGGATATAGC 

: GAAAGGGGAGGAGTGGCAAA 

: CGCACGCGGTGTTAAGAAAT 


G 


Boc o o" oH oo omo oo Ucuo Uu uocuuuuuu'uu'utu'u*u"usuudu"u sumus 


个 SNP 在 抗 性 个 体 和 易 感 个 体 中 等 位 基因 频率 , 筛 
选 到 2 个 候选 SNP 位 点 等 位 基因 频率 在 抗 性 个 体 
与 易 感 个 体 中 存在 差异 显著 性 :C11378366T( 所 在 
染色 体 骨 架 为 NW_003378111. 1) 和 T14570310C 
(所 在 染色 体 骨 架 为 NW. 003378155. 1) ,分别 位 于 
第 2 号 染色 体 的 第 1 108 857 位 和 第 11 号 染色 体 的 
第 156 247 位 。 分 析 测 序 峰 图 ,发现 这 两 个 SNP ЭУ) 
存在 3 种 基因 型 : 纯 合 C/C, TT 及 杂 合 C/T( 图 1)。 
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图 1 意大利 蜜蜂 染色 体 DNA 非 编 码 区 与 幼虫 
抗 白垩 病 相 关 的 SNP 测序 示意 图 


Schematic diagram of sequencing SNPs associated 





Fig. 1 
with larva chalkbrood resistance in the non-coding region 
of DNA of chromosomes of Apis mellifera mellifera 
А: 位 于 第 2 号 染色 体 的 SNP (С11378366Т) (所 在 染色 体 骨架 
为 NW_003378111.1) 的 基因 型 Genotype of the SNP C11378366T 
(scaffold: NW_003378111. 1 ) on chromosome 2; B: 位 于 第 11 号 
染色 体 的 SNP (T14570310C) (所 在 染色 体 骨 架 为 NW_ 
003378155. 1) 的 基因 型 Genotype of the SNP T14570310C 
(scaffold; NW_003378155. 1) on chromosome 11. 顶部 2 图 为 TT 
纯 合 个 体 的 峰 型 图 ,中 间 2 图 为 CT 杂 合 个 体 的 峰 型 图 , 底 边 2 
图 为 CC 纯 合 个 体 的 峰 型 图 。 


peak pattern of TT homozygous individuals, the middle two figures 















































The top two figures represent the 


represent the peak pattern of CT heterozygous individuals, and the 
bottom two figures represent the peak pattern of CC homozygous 


individuals. 


SNP 位 点 C11378366T 在 抗 病 个 体 中 C 等 位 基因 的 
频率 显著 高 于 了 T 等 位 基因 的 频率 ,而 在 易 感 个 体 中 
T 等 位 基因 的 频率 高 于 С 等 位 基因 的 频率 ;SNP 位 
点 T14570310C 在 抗 病 个 体 和 易 感 个 体 中 的 等 位 基 
因 频 率 则 与 SNP. 位 点 C11378366T #5 (192). F 
方 分 析 表 明 : 这 两 个 SNP 在 易 感 病 幼虫 与 抗 病 幼 虫 
中 均 存 在 差异 显著 性 (P<0.01)。 为 了 进一步 验证 





























结果 的 可 信 度 ,我 们 将 前 期 发 现 编码 区 与 白垩 病 抗 
性 相关 的 SNP(C2587245T) 在 意大利 蜜蜂 12 个 抗 
病 个 体 和 12 个 易 感 个 体 中 进行 验证 ,发 现在 抗 病 个 
体 中 C 等 位 基因 的 频率 高 于 T 等 位 基因 的 频率 ,这 
与 我 们 前 期 的 结果 是 一 致 的 。 
2.2 意大利 蜜蜂 染色 体 DNA 非 编 码 区 SNP 位 点 
的 扩大 群体 验证 

为 了 进一步 验证 SNP 的 可 靠 性 ,意大利 蜜蜂 抗 
病 样本 与 易 感 病 样本 各 再 取 90 头 , 用 于 再 次 验证 
SNP 位 点 C11378366T 和 SNP 位 点 T14570310C。 结 
果 发 现 ( 图 3) :SNP 位 点 T14570310C 中 TT 等 位 基因 
在 抗 病 个 体 中 基因 频率 为 0.746269 ,而 T 等 位 基因 
在 易 感 病 个 体 中 的 基因 频率 为 0.564706,T 等 位 基 
因 在 抗 病 个 体 中 的 基因 频率 比 在 易 感 病 个 体 中 高 ， 
且 抗 病 幼 忠 与 易 感 病 幼虫 中 的 基因 型 分 布 仍然 表现 
出 显著 性 差异 (P<0.01) ;而 SNP 位 点 C11378366T 
在 抗 病 幼虫 与 易 感 病 幼虫 中 基因 型 分 布 没有 显著 
性 ,因此 ,该 SNP 被 淘汰 掉 。 同 时 ,我 们 发 现 SNP 位 
点 C2587245T P C 等 位 基因 频率 在 抗 性 个 体 较 高 ， 
而 T 等 位 基因 频率 在 易 感 个 体 较 高 ,这 与 前 期 的 结 
果 也 是 一 致 的 。 
2.3 生物 信息 学 分 析 

由 NCBI 数 据 库 对 SNP 位 点 T14570310C 进行 
定位 查找 ,发 现 该 SNP 位 于 10С100578413 基因 5' 
端的 非 编 码 区 SNP, 



































3 讨论 


本 研究 利用 团队 前 期 筛选 出 与 白垩 病 相 关 的 意 
大 利 蜜蜂 染色 体 DNA 非 编 码 区 SNP ,结合 Holloway 
等 (2013 ) QTL 结果 中 第 2 和 11 号 染色 体 与 蜜蜂 幼 
虫 抗 白垩 病 显 著 相关 ,我 们 利用 PCR 测序 的 方法 对 
第 2 和 11 号 染色 体 DNA 非 编 码 区 55 个 SNP 进行 
验证 , 最 终 筛 选 出 了 一 个 位 于 第 11 号 染色 体 上 
LOC100578413 基因 非 编 码 区 的 SNP 位 点 
T14570310C 与 幼虫 抗 和 白垩 病 性 状 相关 ,该 位 点 在 抗 
病 个 体 中 的 了 等 位 基因 基因 频率 更 高 。 因 此 ,SNP 
位 点 T14570310C 为 T/T 基因 型 的 个 体 可 以 判定 为 
抗 性 个 体 , 而 将 СИС 基因 型 的 个 体 作 为 易 感 个 体 ， 
该 SNP 判定 个 体 抗 病 性 或 易 感性 的 结果 与 前 期 编 
人 码 区 SNP 位 点 С2587245Т 的 结果 是 一 致 的 。 后 期 
可 以 利用 该 SNP 筛选 基因 型 为 TAT 基因 型 的 处 女 
王 , 利 用 CO, 处 理 , 可 使 其 产生 未 受精 卵 , 进 而 发 育 
成 T 基因 型 的 雄 蜂 ,采集 性 成 熟 雄 蜂 的 精液 ,通过 
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SNP C11378366T SNP T14570310C SNP C2587245T 
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Fig. 2 Distributions of C and T allele frequencies of SNPs of Apis mellifera ligustica in preliminary verification 
SNP С2587245Т; 编码 区 与 白垩 病 抗 性 相关 的 SNP (SNPs associated with chalkbrood resistance in coding regions) ; Ро: SNP 位 点 处 C 等 位 基因 频 

















率 Frequency of C allele in the SNP loci; Рр: SNP 位 点 处 了 等 位 基因 的 
上 星 号 示 组 间 差 异 显著 ( 卡 方 检验 , Р <0.05). The asterisk above bars 


抗 病 组 DIRA 抗 病 组 


10 
9 


oo 
о o o 


7 
6 
5 


等 位 基因 频率 (%) 
Frequency of alleles 
ы 
5 со DO С 


Y 
© 


2 


= 
о o 


Wis Frequency of T allele in the SNP loci. 图 3 [F] The same for Fig. 3. Ж 





indicates significant difference between groups ( chi-square test, P «0.05). 


B», EHP 


271 БЕНЕН 


易 感 组 抗 病 组 易 感 组 


T 


Resistant group Susceptible group Resistant group Susceptible group Resistant group Susceptible group 


SNP C11378366T SNP 


Т14570310С SNP С2587245Т 


3 意大利 蜜蜂 扩大 群体 验证 中 SNP 的 C 和 T 等 位 基因 频率 
Fig. 3 С and T allele frequencies of SNPs in the validation of expanded population of Apis mellifera ligustica. 


E ns 和 三 星 号 分 别 表 示 组 间 无 显著 性 差异 (P >0.05) 和 差异 极 显 著 (P <0.001) ( 卡 方 检验 ) ns and triple asterisk above bars indicate no 


significant difference (Р >0.05) and extremely significant difference (Р «0.001) between groups, respectively, by chi-square test. 











人 工 授精 的 方式 ,进而 可 以 培养 大 量 抗 白垩 病 蜂 群 。 
同时 由 于 白垩 病 发 病 的 初期 没有 明显 症状, 一 旦 出 


现 明显 症状 ,很 快 就 会 导致 幼虫 死亡 ,该 分 子 标记 为 
白垩 病 早期 的 诊断 和 培育 白垩 病 抗 性 蜂 种 提供 新 的 
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检测 方法 。 

关于 SNP 位 点 T14570310C, 它 所 在 的 基因 
LOC100578413 在 数据 库 中 没有 注释 ,属于 未 知 功 能 
基因 的 非 编码 区 SNP。 非 编码 DNA 是 指 不 编码 蛋 
白质 的 DNA 序列 ,包括 ncRNA 和 miRNA 等 ,能 
通过 自身 转录 、DNA 甲 基 化 、 组 蛋白 修饰 . 重 塑 染色 
质 等 途径 调控 基因 表达 ,影响 多 种 生物 学 过 程 ( 秦 
丹 和 人 徐 存 挫 ，2013 ) 。miRNA 与 生物 的 生长 发 育 、 
免疫 细胞 凋 亡 、 癌 症 等 都 存在 重要 的 关系 (Shomron 
et al., 2009; Zhang et al., 2013) 。 在 黑 腹 果 蝇 中 ， 
miR-14 可 作用 于 Drice 和 sugarbabe 基因 调控 细胞 凋 
亡 及 物质 代谢 (Xu et al., 2003; Varghese et al., 
2010); Æ X Æ, miR-33, miR-190, miR-276 和 
miR-7 n[ fe Js] fa ЛЕНІ ЕНУ Ж ( Burgler and 
Macdonald, 2005) ; ЕЖЕН, 发 现 采 集 蜂 头 部 的 
miR-2796 可 能 与 蜜蜂 的 神经 发 育 有 关 ( Liu et al., 
2010)。 同 时 非 编码 区 DNA. 还 会 影响 疾病 , 例如 
NBN 基因 3' 端 非 编码 区 1510464867 与 胃癌 风险 相 
关 ( 孙 屏 等 , 2019) ;APC( 结 肠 腺 瘤 性 息肉 病 ) 基因 
3 ' 端 非 编 码 区 rs1804197 多 态 性 与 结 直肠 癌 易 感性 
相关 ( 陈 志 鹏 等 , 2017) ;表皮 生长 因子 5' 端 非 编码 
区 rs444490 位 点 多 态 性 与 吧 母 细胞 瘤 易 感性 及 相 
关 ( 周 洲 和 朱 房 勇 , 2016 ) 。 基 因 非 编码 区 的 SNP 
要 比 编码 区 的 SNP 多 ,可 以 影响 蛋白 质 的 结构 和 功 
能 ,影响 转录 因子 和 顺 式 作用 元 件 ,如 启动 子 、 增 强 
子 、 绝 缘 子 、 沉 默 子 等 的 结合 导致 疾病 的 发 生 
( Shomron et al., 2009; 戎 蛤 等 ,2016)。 如 在 人 类 
医学 方面 的 研究 证 实 beta-globin ( HBB) 基因 增强 子 
的 缺失 和 重组 可 导致 地 中 海 贫血 , Sonic hedgehog 
(SHH) 基因 limb 增强 子 的 点 突变 导致 轴 前 多 指 症 
( Visel et al., 2009) 。 此 外 ,还 有 研究 发 现在 人 类 染 
色 体 上 一 个 长 度 为 136 kb 片段 的 着 丝 粒 和 端 粒 区 
域 存在 一 些 SNP 位 点 与 前 列 腺 癌 和 结肠 癌 相 关 , 部 
分 位 点 位 于 高 保守 序列 的 增强 子 上 ( Alexander. et 
al., 2010)。 尤 其 是 一 些 位 于 基因 非 编 码 区 的 SNP 
还 能 通过 影响 mRNA 的 结构 导致 一 些 疾病 的 易 感 
性 ,如 丙 型 肝炎 (Lu et al., 2015) 。 研 究 表明 ,位 于 
miRNA 靶 基 因 结 合 位 点 的 一 个 或 多 个 SNP 位 点 与 
许多 人 类 疾病 在 流行 病 学 上 显著 相关 。 目 前 研究 表 
明基 因 非 编码 区 SNP, 如 miRNA Му) SNP 与 眼 部 
疾病 的 易 感性 状 有 关 ( 周 洲 和 朱 房 勇 , 2016 ) ІЗ 
此 ,SNP T14570310C 在 白垩 病 抗 性 个 体 和 易 感 个 体 
之 间 存 在 显著 差异 ,暗示 该 SNP 可 能 通过 某 种 途径 
导致 蜜蜂 幼虫 对 病原 的 抵抗 能 力 产 生 影 响 。 
































































































































目前 ,基因 组 编辑 技术 CRISPER-Cas9 已 在 蜜 
蜂 中 应 用 ,为 基因 功能 验证 提供 了 有 力 的 技术 保证 。 
我 们 后 期 可 以 对 SNP. T14570310C 位 点 进行 编辑 ， 
深入 研究 该 SNP 在 白垩 病 抗 性 机 制 中 的 作用 。 
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